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Abstract 



We discuss the construction of basis states for Hamiltonian QCD on 

the lattice, in particular states with dynamical quark pairs. We cal- 

,JZj . culate the matrix elements of the operators in the QCD Hamiltonian 

between these states. Along with the "harmonic oscillator" states in- 
troduced in previous pure SU(3) work, these states form a working 
basis for calculations in full QCD. 



1 Introduction 

As detailed in papers with Bronzan [0, [|, we have been studying lattice 
QCD using an operator and states approach rather than the usual Monte 
Carlo simulation of the Feynman path integral. As we have emphasized, this 
requires the use of a suitable set of single degree-of-freedom (DOF) states 
on the SU(3) manifold. Here "suitable" means that matrix elements of the 
QCD Hamiltonian must be calculable in closed form using these states, and 
that the states have tunable parameters so they can model the QCD wave 
functions at all values of the coupling constant. 

We have introduced a "harmonic oscillator" basis of states as one which 
satisfies the above criteria and is therefore useful in making calculations in 
Hamiltonian QCD. Since these states describe only gauge degrees of freedom, 
they are actually designed for studying only a pure SU(3) gauge theory. 

In this paper we will introduce a similarly suitable set of dynamical- 
quark basis states to complement the harmonic oscillator states. We will 
discuss how we construct a quark ground state at arbritrary values of the 
coupling constant, create quark excitations, and calculate matrix elements 
of the QCD Hamiltonian between these states. Together with the harmonic 
oscillator states, then, we have a complete basis for the study of full QCD. 

The layout of this paper is as follows. In Section 2 we briefly review 
our work on the construction of harmonic oscillator states. In Section 3 we 
describe the construction of the quark vacuum as well as the states which 
include virtual quark pairs. In Section 4 we discuss the combination of the 
harmonic oscillator states with the quark states to form a basis for simulating 
full QCD. In Section 5 we derive the matrix elements of the operators of the 
QCD Hamiltonian between these states. We conclude in Section 6 with a 
discussion of the usefulness of these states. 



2 Review of pure gauge states 

Each gauge degree of freedom (link variable on the lattice) can be described 
by a wave function which is a "harmonic oscillator" state on the SU(3) man- 
ifold. The harmonic oscillator states are derived from a "Gaussian" state on 



the manifold. This Gaussian can be written 

^ ai (a } t) = Y,d(p,q)e- tx ^ X ^\ (2.1) 



p,q 



where d(p,q) — (p + l)(q + l)(p + q + 2)/2 is the dimension of the (p, q) 
representation of SU(3), X(p, q) = (p + g) 2 /4 + Qo — g) 2 /12 + (p + q) is the 
quadratic Casimir eigenvalue, and y^'^ is the character, a = (a 1 , a 2 , . . . , a 8 ) 
is a parameterization of the adjoint representation of SU(3). The parameter 
t controls the width of the Gaussian, which is centered at a±. 

To make connection with a more familiar object, note that the corre- 
sponding wave function on the flat, three-dimensional manifold R 3 would 
be 

e -(x- Xl ) 2 /4t 

where x and xi are ordinary vectors. On this flat manifold, harmonic os- 
cillator states can be generated by applying a polynomial in the operators 
—i-Q^— to the wave function, then setting x x to 0. On the SU(3) manifold, the 
corresponding operators are J^Lai&i), the generators of SU(3) in differential 
form.[] 

The harmonic oscillator states that we have formed on the manifold can 
be labeled by the number of SU(3) indices on the state. Thus, the zero-, 
one-, and two-color states are 
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(2.3) 



(2.4) 

1 Whcn applied to a representation of SU(3), Jh a {. a \) has the effect of left-multiplying 
the representation by the a th generator. There are also operators jR a {a\) which right- 
multiply the representation and could be used to generate a slightly different set of har- 
monic oscillator states. 



In the papers we have cited we described the orthogonalization of these states 
as well as the calculation of matrix elements of the QCD Hamiltonian between 
them. We will not need those details here, however. 

These single-DOF wave functions must now be pieced together to form 
multi-DOF states suitable for computations in Hamiltonian QCD. Each multi- 
DOF state will be a product of single-DOF states. The most basic state is 
one with no excitations. Specifically, each degree of freedom is a zero-color 
state: 

|#>= |0000... > . (2.5) 

The next simplest states are those with just one excitation: 

|^>=|0 a 000--->. (2.6) 

However, these states are not SU(3) singlets, as evidenced by the color index 
which is left hanging. In order to maintain global SU(3) invariance, we must 
sum over all indices. Therefore, the next possible states are 

|^>=<U|0 a b 00--->, (2.7) 

where 5 a b is a Kronecker ^-function. There are of course similar states where 
the excitations lie on different degrees of freedom. We call these one-pair 
states. In analogous fashion, we can form two-pair states, three-pair states, 
and so on. 

Arbitrarily intricate states can be constructed by using highly colored, 
single-DOF states along with complicated coupling schemes. For instance, 
the following is a valid (that is, gauge-invariant) state: 

|^>=<UTr(A c A d A e )|0 ac bd 00 e --->, (2.8) 

where the Aj's are Gell-Mann matrices. 

The states described above provide the basis for a pure SU(3) gauge 
theory calculation, and they will likewise provide the basis for the gauge 
sector of our full QCD calculation. In the next sections we will describe how 
to form similar states that include dynamical quarks, to complete our basis. 



3 Quark sector states 

To formulate our quark states, we use staggered fermions ||. (We follow 
closely the notation of Banks et al. PJ). Specifically, the quarks are created 



and destroyed by applications of single-component fields xl( s ) an d X«( s )> 
where the index % signifies color in the fundamental representation of SU(3), 
and the argument s = (x, y, z) specifies a site on the lattice. The operators 
act on a simple Fock space, and any quark configuration can be specified by 
giving the location of the quarks on the lattice. 

3.1 Strong-coupling quark vacuum 

Just as we did for the gauge states, we would like to use quark states which are 
flexible enough to treat both the strong- and weak-coupling regimes. Toward 
this end, we will begin with a description of the strong-coupling vacuum, 
then discuss how to obtain the weak-coupling vacuum from it. We will then 
describe the excitations of the vacuum. 

At very strong coupling the vacuum state is the "checkerboard" state, 
in which alternating sites are fully occupied |j]. Two of these states can 
be formed, depending on which sublattice is filled (the "red" squares or the 
"black" squares of the checkerboard). The two states are 

\ll>even> = *[[ ^CijkXi( a )Xj(»)Xk( S ) \E >, 

even sites 

\^odd> = II g e yfeXi( s )xJ(s)xl(s) \E>, (3.1) 

odd sites 

where \E> is the empty lattice, and even/odd means that x + y + z = ±1.0 
These two states transform into one another under a chirality transforma- 
tion. We can form two degenerate vacua, the symmetric and antisymmetric 
combinations of these two states, which are eigenstates of the chirality oper- 
ator. We have shown numerically that as we decrease the coupling constant, 
the true vacuum is the symmetric combination, so we use 

\0>strong= \^even> + \^odd> (3.2) 

as the vacuum state at strong coupling. 



2 One must be careful to maintain a fixed ordering for the product over the sites, lest 
minus sign errors creep in. Although we will not display it explicitly, we assume that the 
sites have been ordered in some appropriate manner. 



3.2 Weak-coupling quark vacuum 

We derive the weak-coupling vacuum from the strong-coupling vacuum by 
means of a projection operator. We write 



\0>weak= lim e aHw \0>strong ■ (3.3) 

H w is the weak-coupling limit of the quark Hamiltonian, 

^ = ^Ex t (s)M ss , X (s'), (3.4) 

s,s' 

where 

M ss > = X^s,s'+n + <5s,s'-nMn), (3.5) 

T](x) = (— l) z , r)(y) = (— l) x , r)(z) = (— l) y , and a is the lattice spacing. 
(There is a suppressed color index in this expression.) Note that in weak 
coupling the quarks decouple from the gauge degrees of freedom, and H w is 
equal to the free-field quark Hamiltonian. 

The projection works as follows, a is an adjustable parameter which 
governs the strength of the projection. When a = 0, there is no projection, 
and we still have the strong-coupling vacuum. As a becomes large, the 
projection operator damps out the pieces of |0 > strong which have high energy 
in the weak-coupling limit, leaving intact the lowest energy weak-coupling 
state. So, for large a, this operator does exactly what we want: it projects 
out the weak-coupling vacuum.^ In practice, a = 1 will give nearly complete 
projection. As we vary a from to 1, we obtain a vacuum state which 
interpolates smoothly between the g — > oo and g = vacua. We therefore 
hope to obtain a good approximate quark vacuum for arbitrary values of the 
coupling constant by using a as a variational parameter in our simulations. 
We will designate this approximate vacuum |0 9 >. 

3.3 Quark excitations 

Having established a quark-sector vacuum, we must now understand how to 
create excited quark states. This is much easier and more familiar than it 



3 Of course, we are assuming here that the strong- and weak-coupling vacua are not 
orthogonal to each other. This can be demonstrated a posteriori. 



was for the gluon states, where we were dealing with the complicated SU(3) 
manifold. Here, we have a simple Fock space, and we merely need to decide 
which states from that Fock space to include in our basis. 

We can exclude a large number of states by observing that the "checker- 
board" vacuum state is half filled and that the Hamiltonian preserves quark 
number. We are interested in studying the glueball spectrum, which has the 
same quantum numbers as the vacuum, so there is no reason to include any 
quark state that isn't half-filled. If we wish to study the spectrum of particles 
with other quantum numbers, we likewise limit ourselves to those sectors. 

We can obtain a new half-filled state directly from the vacuum state by 
applying one annihilation operator and one creation operator: 

|*>=xl(s) Xi (s')|O g >. (3.6) 

Recall that the index i indicates color, in the fundamental representation 
of SU(3), reflecting the correct group properties of the quarks. Notice that 
we are summing over i, thereby maintaining global gauge invariance in the 
quark sector in an analogous way to the gauge sector. 

The state shown above corresponds to exciting one virtual quark pair 
out of the vacuum. Clearly, we can form states with an arbitrary number of 
virtual pairs by applying the appropriate number of creation and annihilation 
operators. One convenient feature of this scheme is that we can control how 
many such pairs we wish to have in the simulation simply by choosing which 
basis states to include. We will comment on this in the final section. 



4 Combining the gauge and quark states 

Sections 2 and 3 described our construction of states for the gauge and quark 
degrees of freedom in our simulations. The simplest way to form a combined 
state suitable for a full QCD simulation is to make a direct product of one 
state from the gauge sector with one from the quark sector. A typical state, 
then could be written 

|*>= SabScd\<f>a<f>b<f>c<f>d<f>---> ® xl(Sl)Xi(s / 1 )x j t (s 2 )Xi(s / 2 )|0 g > . (4.1) 

Let us emphasize again that all of the SU(3) indices are contracted, ensuring 
that global gauge invariance is maintained. Note in particular that the indices 

6 



are summed separately within the two sectors, indicating that the gauge and 
quark sectors are separately invariant under the global SU(3) transformation. 
There is a second way to form combined states, in which the two sec- 
tors are not separately invariant. To do this, we need an object which can 
intertwine the SU(3) adjoint representation of the gauge sector with the fun- 
damental representation of the quark sector. The simplest way to do this is 
with the Gell-Mann matrices, which carry both three- and eight-dimensional 
indices. A simple state of this form is 

|M>>= |0 a 0000> <g> xi(s)A^(s')|O g > . (4.2) 

This completes the description of our basis states for full QCD calcula- 
tions on the lattice. In the next section we will show how to compute the 
matrix elements necessary for our simulations. 



5 Matrix elements 

The matrix elements of our combined states are sums of products of gauge 
matrix elements with quark matrix elements. The results of our gauge sector 
calculations have already appeared in ||, and we will not repeat them here. 
The most basic quark matrix element that we must calculate is the overlap 
of our projected vacuum, <0 g |0 g >. This, in turn, can be broken down into 
sums of matrix elements of the projected checkerboard states: 

<0 g |0 g > = <i)even\e~ 2aHw \i)even> + <^even\e~ 2aHw \i)odd> 

+ < *Podd\e- 2aHw |Vw > + < ^odd\e- 2aHw \iJodd > ■ (5.1) 

We will explicitly calculate only the first of these four nearly identical ex- 
pressions. Recall that 



|VWn> = II Q e ijkxl(s)x](s)xi(s)\E> 

even sites 

= \x\MxlMxlM\- ■ ■lx\(sM)xi(sM)xi(sM)]\E>, (5.2) 



where M = ^- is the number of even sites on our N x N x N lattice,^ 
and Si, . . . , Sm specify the locations of those sites. First, we segregate the 



We assume N is an even number. 



creation operators by color. This involves an even number of permutations, 
so no minus sign is picked up. 

\^even>= [Xl(Sl) • • • Xl (%)] \xU S l) ' * * X^M )] [Xs( S l) * ' -x1^m)]\E> • 

(5.3) 
Likewise, we can segregate H w into a product of single-color factors: 

= ^E[x t l( s ) M ss'Xl(s / )^-xl(s)M ss , X2 (s')^-xks)M ss , X 3(s , )]■ (5.4) 

s,s' 

Since the separate color factors commute, we can write 

e -2aH w _ e -2aHl e -2aHl e -2aHl (5.5) 

and thereby write our matrix element in color-factorized form: 

<^even\e- 2aHw \^ eV en> = < E\ [ X l (s M ) " ' " Xl (Si)]e" 2aH i [ X \ (Si) • • • XlM] 

x [ X2 (sm) • • •X2(si)]e- 2aff ™[ X 5(si) • • "X 2 (sm)] 
x [x 3 (sm) • • •X3(si)]e- 2 ^[ X3 (si) • ■■xl(s M )]\E> . 

(5.6) 

We now want to commute the exponential factors all the way to the left, 
where we will be able to take advantage of the fact that 

<E\e~ 2aH ™ =<E\. (5.7) 

To do this, we will need the relation 

X(s)e-**- = ]>> s , s ,(<T)e-^ X (s'), (5.8) 

s' 

where 

K S:S ,(a) = $>(s;Me CTA ( k <*V(s';M (5.9) 

and the functions </>(s;k, i) are the eigenf unctions of M ss /. We prove this 
relation in the Appendix. 



Commuting each of the exponential factors to the left gives 



M 



< ileven\e- 2aH ™\i>even> = < E\ J{[£K SiA {2a)xM)}x\{^) ' ' ' xlM 

=1 s' 

M 

IIE^-i,<( 2 «)X2(s;)]xS(si)---xJ(sM) 



i=\ 
M 
X 

i=l s ; 

M 
X 

i=l s 



IIE^-i,<( 2 «)Xs(sJ)]xJ(si)"-xS(sM)|S>. 



(5.10) 



The M sites s' x , . . . , s' M can be chosen in any way from the even sites, but 
unless they are a permutation of Si, . . . , Sm, the contribution to the matrix 
element is zero. Thus, we get 

<1peven\e~ 2aHvJ \^even> = <E\ J^ K si,s P(1) ' ' ' #s M) s P(M) ep(l)...p(Af) 

perms P 

x [Xi(si) • • • Xi(sm)] [Xi(si) • • • Xi(sm)] 
x {Similar factors for other colors} \E> 
= I Z^ E sl!Sp(1) ■ • • K SMjSp(M) \ . (5.11) 

perms P 

Introduce the M x M matrix L EE (a), whose (m, n) th element is given by 

LH(a) = K Sm , Sn (2a). (5.12) 

The EE refers to the fact that we are calculating the even-even matrix ele- 
ment. We can then write our even-even matrix element in its final, decep- 
tively compact form: 

< i>even\e- 2aHw \ll>even>= [det L EE {a)f . (5.13) 

The even-odd, odd-even, and odd-odd pieces of <0 (? |0 (? > can be calculated 
in analogous fashion, the only difference being the composition of the corre- 
sponding L matrix. We therefore have completed the calculation of the norm 
of the quark vacuum: 

<0 q \0 q >= [detL EE (a)] 3 + [det L EO (a)] 3 + [det L OE (a)] 3 + [det L°°(a)] 3 . 

(5.14) 

9 



Now that the machinery is set up, the remaining matrix elements require 
much less effort. The next simplest matrix element to calculate is 

<0,|xi(so)xJ(s[,)|O ff >. (5.15) 

We need not consider operators XX or X X > which change quark number and 
therefore have trivially vanishing matrix elements. (As we have mentioned 
already, the QCD Hamiltonian preserves quark number, and therefore does 
not include operators of this type.) Similarly, if % ^ j, then the matrix 
element vanishes because there would be a mismatch in the number of quarks 
of a given color. 

The calculation of this matrix element proceeds exactly as the overlap 
calculation, except that one of the colors (suppose it's color 1) will have the 
extra operators Xi( s o) and Xi( s d) to commute through the exponential. This 
gives color 1 a contribution of 

J2 K s s (-a)K Ks , o (-a) <£|[xi(s M )---Xi(si)xi(s )] 

So, So 

x e- 2 ^[x t i(s , )x t i(si)---X t 1 (s M )]|^> -(5.16) 
If we define the extended (M + 1) x (M + 1) matrix 



L (s ,s a 



( ^.sf, 


^s ,si(a) • 


' K S(hSM (a) \ 


^si,s (a) 


^s 1)Sl (2a) • 


■ K SuSM (2a) 


V K SMK (a) 


K SM>Sl (2a) • 


■ K SM>SM (2a) ) 



(5.17) 



then we can again write a compact expression for our matrix element, 

<0,|Xi(so)x5(so)|0,> = (%detL EE (s ,s(,)[detL EE (a)] 2 

+ EO + OE + 00 contributions. (5.18) 

The factors of detL EE (o;) are the contribution of the two unmodified colors. 
Matrix elements involving more XX^ pairs can be calculated in analogy 
with the above derivation. Each additional pair requires the introduction of 
a further extended matrix, for the case in which all the x's act on the same 
color. For instance, the next extension of the matrix comes from 

<0 (? |xi(s )xi(soo)x t i(so)x t i(s[ )0 )|0 (/ > = detL EE (s ,s[ ) ;soo,s[ )0 )[detL EE (a)] 2 

+ EO + OE + 00 contributions, 

(5.19) 



10 



where 



( "s 00 ,s£, 
"s o,sf, 




^80,81 (a) • 


•• ^s o,s Af («) \ 
■■ ^s ,sm(") 


^ 11S ' 00 («) 


K S!,S' (01) 


K SltSl {2a) ■ 


■■ ^ SlJ s M (2a) 


V ^s Mj s' 00 (a) 


K SM ,s' ( a ) 


^s M ,si(2a) • 


•• #s M ,s M (2a) ) 



L ( s o, s oi s oo, s oo) 



(5.20) 

We have calculated matrix elements of this type for up to four xx^ pairs, 
which is the requirement for simulations with two dynamical quark pairs.0 



6 Summary 

We have introduced a set of basis states that is suitable for making calcula- 
tions in an "operator and states" approach to the study of lattice QCD. The 
previous section detailed the computation of matrix elements of the quark- 
sector states. In previous work we calculated matrix elements for pure gauge 
simulations. Taken together, these matrix elements provide us the tools to 
construct a Hamiltonian matrix on a basis of states in full QCD. The eigen- 
values of this matrix are then estimates of masses in the QCD spectrum. 

We would like to emphasize the particular power that this method has 
with respect to testing the quenched approximation. Unlike Monte Carlo 
simulations of the Feynman path integral, dynamical quarks pose no partic- 
ular problems to our method. The states which include quark pairs are "just 
another configuration." Moreover, by choosing which states to include, we 
can limit the simulation to having as many virtual quark pairs as we wish. 
By comparing runs with and without quark pairs, we hope to be able to 
provide some insight into the quenched approximation. 

We would like to acknowledge the contributions made to this work by 
J. B. Bronzan. 



5 If one fixes the gauge by our current method |5|, the QCD Hamiltonian can contribute 
up to two pairs of operators, rather than just one. In this case we would need two more 
"extensions" of the L matrix to do the equivalent simulations. 
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A Appendix 

To prove the relation 

X (s)e-^ = £*M'(*)e~**.x(s / ), (A . 1} 

s' 

we first need to diagonalize H w . Recall that H w is bilinear in the quark 
creation and annihilation operators: 

ff M = ^ X f (#Ws'), (A.2) 

s,s' 

where 

M ss > = X^s.s'+n + ^s,s'-n]^(n). (A. 3) 

fi 

Diagonalizing M ss / amounts to solving the Dirac equation on the lattice. We 
introduce quark mode operators, which annihilate "plane-wave" states: 

X (M = $>(s)<Ks;M. (A.4) 

s 

The 0(s; k, i) are the eigenf unctions of M ss r. 

£ M ss ,0(s'; k, i) = A(k, i)0(s; k, z), (A.5) 

s' 

where k is the lattice momentum, i enumerates the modes, and A(k, i) is the 
eigenvalue of the mode.0 On a lattice with NxNxN sites, 1 < k x , k y , k z < y 
and 1 < % < N 3 . The derivation of the eigenfunctions is straightforward but 
rather tedious. We have carried it out explicitly only for a 2 x 2 x 2 lattice. 
The weak-coupling quark Hamiltonian is (by design) diagonal using these 
new operators: 

s,s' 

= ^EEE#. *)<K s ; k, *)m ss ,0(s' ; k', i') x (k f , o 

s,s' k,k' i,i' 

= l^A(k,z) X t( k ,i) x (k,i). (A.6) 

Za k,i 



D The plane-wave operators also have a color index, which we suppress here. 
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Next, we define the expression 

S(a) = e aH ™x(Ki)e~ aHw . (A.7) 

Then 

--e ffH -[H w , X (k,i)}e- ffH -. (A.8) 



da 
Substituting the diagonal form of H w we just derived, we have 

[H w ,x(k,i)] = -\(k,i) x (k,i). (A.9) 

Therefore, 

^M = -\(k, i )S(a), (A.10) 

which implies that 

S(a) = constant x e ^ A(M . (A. 11) 

Noting that 5(0) = x(k, i), we have 

S(a) = e-° x ^ X (Kt), (A.12) 

giving us the relation 

e aHw X {K i)e~ aHw = e- <TA(k>i) x(k, i). (A.13) 

If we now multiply both sides of this by 0(s;k, i), sum over k and i, and 
change back to configuration space variables x( s )> we obtain 

e^ X (s)e-^ = £ £ 0(s; k, *) e -^0(s'; k, i) x (s'). (A.14) 

s' k,i 

Substituting the variable fT, and left-multiplying by e^ crHw , we get the rela- 
tion we set out to prove, 

X (s)e- ffH ™ = £* M '(*)er" ir -x(s / )- (A.15) 



s' 
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